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We present a simple and intuitive picture for the deconfinement of quarks and gluons at finite 
temperature: as the temperature increases, QCD behaves like QED at T = 0. We show this by 
calculating the QCD coupling constant as a function of the temperature and of the external momenta 
used to probe quarks and gluons. 



The study of QCD at finite temperature and density is of crucial importance for the understanding of hadron 
formation in the early universe: as the temperature decreased, quarks and gluons must have experimented a phase 
transition from a quark-gluon dominated universe to a hadron dominated universe. Efforts to recreate the quark <-> 
hadron transition are currently being fully implemented at RHIC, where it is hoped that clear signals of a quark-gluon 
■ plasma phase will be detected. 

It is also currently believed, thanks to lattice simulations of finite temperature QCD [Q, that a deconfined phase is 
i—i ■ expected to exist at temperatures above (and around) 0.2 GeV. However, we still do not have a clear picture of how 
I ^ | QCD itself changes from a confining to a non-confining theory. That is, the quarks which were originally enslaved 
inside a region of ~ 1 fm are now allowed to propagate over a larger portion of space, outside the region occupied 
by the original nucleoli. This is the strict sense that the word deconfinement is used througout this work. On the 
other hand, the difference between a non-confining theory, QED in the case, and a confining theory, QCD, resides 
exactly in the fact that the later has a gauge self coupling^. Physically, it may help to think that the deconfinement 
t^J- \ proccess considered here involves a transition from a non Abelian to an Abelian theory. In this sense, we think that 
f**«> ■ in a deconfined phase, QCD should be "less" non Abelian than in a confined phase. 

It is well known that at T — 0, besides a color factor, the QED and QCD f3 functions would be the same if only 
fermion loops existed. The difference in sign between the QED and QCD (3 functions comes, essentially, from the 
gluon and ghost loops in the gluon self energy (plus, of course, the vertex corrections). Alternatively, we notice that 
the two beta functions are the same in the limit that the number of colors, N c , goes to zero. At finite temperature, 
it would be natural to expect that the same kind of limit exists. We will see in this letter that, in fact, QCD at high 
temperature and small external momenta squared behaves like QED at T = 0, a conclusion naturally drawn from the 
QCD (3 function at finite temperature. However, the calculation of this quantity has been plagued with dificulties, 



o 
o 
o 



i 



most of them related to the lack of gauge invariance 1 3|-]i| and the dependence of the renormalized coupling at finite 



D temperature on the vertex chosen to renormalize it (6-^ 

' To ensure the gauge invariance of the (3 function means to ensure the gauge invariance of the vertex and self energy 

^ , functions. Our freedom to choose among the different vertices (quark-gluon, gluon-gluon, ghost-gluon) should not 

■ affect the calculation of the running coupling. In general, we have: 

a s = Z~ l a B , (1) 

_ _ i 

where Z a is the renormalization constant of the coupling. If the gluon-ghost vertex is used for the renormalization, 
then: 

= — (2) 

with Zti, -Zt3, and Zt3ym, the gluon-ghost vertex, the ghost propagator, and the gluon propagator renormalization 
constants, respectively, calculated at finite temperature. Alternatively, we could use the triple gluon vertex: 

Z a = ^. (3) 

T3YM 



1 QED at strong coupling can also have a confining phase However, in our case, we will be considering only the region of 
small external momentum where, as usual, a® ED ~ 1/137. 
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Here, Ztiym is the renormalization constant of the triple gluon vertex. A similar expression exists for the quark 
vertex. 

In a consistent calculation, Z a is independent of the vertex used - at zero temperature this is a well known 
consequence of the Slavnov- Taylor identities. At finite temperature the situation is the same, as long as we restrain 
ourselves to the dominant, gauge invariant, terms as given by the Hard Thermal Loop (HTL) resummation program 
of Braeten and Pisarsky |J. In fact, within this program one can show that the dominant terms obey the Abelian 
Ward identities §,0: 

fc"n M „(fc)=0, (4) 



k fl T fJlUp {k,p,q) = Iiup{q) 



H,p(p). 



(5) 



The HTL approximation consists of disregarding the external momenta in the numerators of the loop integrals, as 
the main contribution to these integrals comes from the region where the internal momenta are of order of T, with T 
taken to be large. Hence, if we compare the gluon, quark and ghost self energies, we see that the later is subleading 
to the formers because it does not have enough powers of internal momenta in the numerator to produce the leading 
(T 2 ) behavior in the temperature. From Eq. (||) we see that the ghost vertex will be subleading with respect to the 
quark and gluon vertices as well. Taking this into account, it follows from Eqs. (||) and (||) that the thermal parts of 
the triple gluon vertex renormalization constant and of the gluon self energy renormalizion constant should be equal: 



{leading T 2 
J T\YM 



{leading T 2 
J T3YM 



(6) 



This relation is valid only for the dominant, gauge invariant, contribution. Z a could also be calculated through the 
quark vertex. In this case there would appear an other constraint for the leading T 2 terms. Instead of Eq. ra), we 
would have that the renormalization constants for the quark vertex and the quark self energy should be the same. 
A direct computation of the transverse part of the gluon self energy gives [[ll 12 1: 
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where k is the four momentum of the external gluon. Only the dominant term in the high temperature expansion 
was written because of its gauge invariance. One can see that Eq. (Q) vanishs for kg << \k\, and that it is reduced to 
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for ko >> \k\, where a s = <? 2 /47r and /it is a mass scale (fij, = \k\ko). Including the T — part, we can calculate the 
renormalization constant for n^, which will depend on T/fi: 
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where /i is the renormalization scale, a is the gauge parameter, and e comes from the dimensional regularization of the 
T = part of the gluon self energy |l3| . The a dependence in Eq. ([To]) is cancelled by the remaining renormalization 
constants in Eqs. ^) and (||). Finally, we notice that the calculation of the triple gluon vertex gives a T 2 //i^ 
dependence for the dominant term in the same form as Eq. (|^) , in a way that we can define a thermal renormalization 
constant which always satisfies Eq. (||). 

The calculation of the thermal j3 function (including the T — and the T ^ parts), for a fixed temperature but 
an arbitrary renormalization point is now straightforward. Using Eqs. ([!]), (|^), and (fl0|), we have: 
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to order a 2 , where B{T 2 , fj, 2 ) = 2/3 ir 2 T 2 /^ 2 . This is the only Renormalization Group Equation for a s because there 
is only one renormalization scale for both the T = and T ^ parts of the renormalized gluon self energy. The 
solution of Eq. (O) is: 
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(12) 



It is helpful to rewrite Eq. (12) in the same format of the T = theory. To this end, we define an effective number 
of colors and an effective number of flavors: 



12 B(T 2 ,Q 2 )-B{T 2 ,Q 2 ) 
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Nr. 
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,B(T 2 , Q 2 ) - B(T 2 ,Ql) 
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With the help of Eqs. ( |l3| ) and ((T^), the expression for the running coupling is written as: 
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The interesting point of this set of equations is that for fixed Qq and T 2 , N^tf < N c and n% > rif as Q 2 — > 0. It 
implies that at some small value of the external momenta squared, the coupling between quarks and gluons becomes 
that of an Abelian theory. To quantify this assertion, we make an explicit calculation of the coupling as a function 
of Q 2 for some fixed values of T. For a s (Qo)j we use the experimental value measured at mz ~ 91 GeV and at zero 
temperature That is, we assume that at such high values of the virtuality of the probe, temperatures of the 

order of 1 GeV are not relevant, an approximation which amounts to set B(T 2 = 1 GeV 2 , Qq w 83 x 10 3 GeV 2 ) w 0. 

In Fig. we show the behaviour of a s for 3 values of the temperature. At T — GeV we have, as usual, that the 
coupling grows rapidly for Q 2 < 1 GeV 2 . However, at T = 0.5 and 1 GeV, a s (Q 2 ) starts to change its behaviour in 
the region around Q 2 — 10 — 20 GeV 2 . Instead of the rapid growth observed at the T = case, for finite T there is 
first an almost Q 2 independence of the coupling, and then it decreses with Q 2 . This behaviour of a s (Q 2 ), for finite T 
and small Q 2 , is that typical of an Abelian theory. In fact, if we define an effective (3 funtion at one loop, 



fit = Hjvj// - \n 



(16) 



we notice that it changes sign at some small value of Q 2 . Figure || tells us that the effective f3 function changes sign 
around 0.05 and 0.2 GeV 2 for T = 0.5 and 1 GeV, respectively. It means that at those values, the theory changes 
from a non Abelian one (P e ^ > 0) to an Abelian one (P e ^ < 0), implying that at high temperature quarks and 
gluons behave like electrons and photons. 

Because the coupling does not grow at small Q 2 and high T , we are allowed to use perturbation theory in this 
region. A direct consequence of this fact is that we can, in principle, fix a s (Q 2 — > 0, high T) using elastic scattering 
of a quark-gluon plasma by an electron beam. 

The results presented here are easily extended to the case of a finite chemical potencial, /i cp . One just has to 
replace T 2 by T 2 + 3fi 2 p / it 2 in the quark loops in the gluon self energy. The qualitative behaviour of a s and j3 e ^ are 
unchanged by the introduction of a chemical potential: the quark and gluon deconfinement at finite T and n cp still 
proceeds through a transition from a non Abelian to an Abelian theory. 
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FIG. 1. The strong coupling constant as a function of Q 2 calculated for 3 differerent values of the temperature. 
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FIG. 2. The effective beta function as a function of Q 2 for two values of the temperature. For large Q 2 , (i e '' tends to 
which is the value of the j3 function at T = for N c = nf = 0. 



